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Abstract 

At large Nc the nucleon can be viewed as a soliton of the effective chiral lagrangian. 
This picture of nucleons allows a consistent nonperturbative calculation of the leading- 
twist parton distributions at a low normalization point. We derive general formulae 
for the polarized and unpolarized distributions (singlet and non-singlet) in the chiral 
quark-soliton model. The consistency of our approach is demonstrated by checking the 
baryon number, isospin and total momentum sum rules, as well as the Bjorken sum 
rule. We present numerical estimates of the quark and antiquark distributions and 
find reasonable agreement with parametrizations of the data at a low normalization 
point. In particular, we obtain a sizeable fraction of antiquarks, in agreement with the 
phenomenological analysis. 
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1 Introduction 



The distributions of quarks, antiquarks and gluons in nucleons, as measured in the 
inclusive deep inelastic scattering of leptons, provides us probably with the largest 
portion of quantitative information about strong interactions. Until now only the 
evolution of the structure functions from a high value of q'^ to even higher values, has 
been successfully compared with the data. It is the field of the perturbative QCD, 
and its success has been, historically, essential in establishing the validity of the QCD 
itself. Unfortunately, the initial conditions for that evolution, namely the leading- 
twist distributions at a relatively low normalization point, belong to the field of the 
nonperturbative QCD, and the success here is still rather modest. 

In this paper we attempt to calculate parton distributions at a low normalization 
point in the limit of large number of colours, Nc ^ oo. Even though in reality Nc = 3, 
the academic limit of large Nc is known to be a useful guideline. At large Nc the 
nucleon is heavy and can be viewed as a classical soliton of the pion field An 
example of the dynamical realization of this idea is given by the Skyrme model Q. 
However, the Skyrme model is based on an unrealistic effective chiral lagrangian. A 
far more realistic effective chiral lagrangian is given by the functional integral over 
quarks in the background pion field M, H : 



exp(iS'eff[vr(x)]) = J D^^Dijexp i^i J (fxil}{i^ - MU^^)^) 

U = exp {iTr''{x)T^) , U^' = exp (i7r'^(x)r"75) = ^-^^U + ^^t^^- (1-1) 

Here ip is the quark field, M is the effective quark mass which is due to the spontaneous 
breakdown of chiral symmetry (generally speaking, it is momentum-dependent) and 
U is the SU{2) chiral pion field. The effective chiral action given by eq. ( |1.1[ ) is 
known to contain automatically the Wess-Zumino term and the four-derivative Gasser- 
Leutwyler terms, with correct coefficients. Therefore, at least the first four terms 
of the gradient expansion of the effective chiral lagrangian are correctly reproduced 
by eq. (|1.1| ), and chiral symmetry arguments do not leave much freedom to further 



modifications. Eq. (1.1) has been derived from the instanton model of the QCD vacuum 
[0, m, which provides a natural mechanism of chiral symmetry breaking and enables 
one to express the dynamical mass M and the ultraviolet cutoff A intrinsic in eq. ( p..l| ) 



through the Aqcd parameter. The effective chiral theory (1^) is valid for the values 
of the quark momenta up to the ultraviolet cutoff A. Therefore, in using eq. (1.1) we 
imply that we are computing the parton distributions at the normalization point about 
A « 600 MeV. It should be mentioned that eq. ( |1.1|) is of a general nature: one need 
not believe in instantons and still use eq. ( |1.1| ). 

An immediate implication of this effective chiral theory is the quark-soliton model 
for baryons of ref. Q, which is in the spirit of the earlier works ||7|, |8[ but without 
the vacuum instability paradox noticed there. According to the model nucleons can 
be viewed as Nc (=3) "valence" quarks bound by a self-consistent hedgehog-like pion 
field (the "soliton" ) whose energy coincides in fact with the aggregate energy of quarks 
from the negative-energy Dirac continuum. Similarly to the Skyrme model large Nc 
are needed as an algebraic parameter to justify the use of the mean- field approxima- 
tion (like one needs large Z to justify the Thomas-Fermi atom), however the 1/Nc 



corrections can be and in some cases are computed [y, IC]. The quark-soliton model of 
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nucleons developed in ref. Q includes a collective-quantization procedure to deal with 
the rotational excitations of the quark-pion soliton. (The quantization of the otherwise 
static solution is necessary to obtain physical baryon states with definite quantum num- 
bers). It enables one to calculate the N and A properties, such as formfactors, A — 
splitting, magnetic moments, axial constants, etc. For a review of baryon properties 
obtained from the model see |11| and references therein. Until now the model lacks 
explicit confinement (though probably it can be implemented along the lines discussed 
in ref. [^]), but it seems to be not so important for the ground state nucleon. 

Turning to the calculation of the nucleon structure functions we note that the 
model possesses all features needed for a successful description of the nucleon parton 
structure: it is an essentially quantum field-theoretical relativistic model with explicit 
quark degrees of freedom, which allows an unambiguous identification of quark as well 
as antiquark distributions in the nucleon. This should be contrasted to the Skyrme 
model where it is not too clear how to define quark and antiquark distributions. The 
advantage of the model can be also seen if one compares it to any variant of the bag 
model. Unfortunately, the bag surface is not described consistently in terms of fields. 
Ignoring the "structure function of the surface" the bag models run into a violation of 
general theorems requiring the complete account for all constituents (see section 8). It 
should be added that all modern fits to the data tend to include antiquarks and gluons 
at a low normalization point, below 1 GeV^ [|^, 14, 15|. 

In this paper we develop the framework for calculating polarized and unpolarized 
quark and antiquark distributions of the nucleon as described by the effective chiral 
theory We check the validity of general theorems, like the sum rules for the 

baryon number, isospin and the total momentum carried by partons, as well as the 
Bjorken sum rule for the polarized distributions. We also derive the expression for the 
Gottfried sum rule and show that its r.h.s. is generally non-zero. 

We show that, from the viewpoint of large Nc, all distributions can be divided 
into "large" and "small" ones. We estimate numerically the "large" distributions: 
the singlet unpolarized distribution and the isovector polarized one, for quarks and 
antiquarks separately Q. The obtained distributions should be, in principle, used as 
initial conditions for the standard perturbative evolution to higher values of where 
one can compare them with the available data. Actually, in this paper we compare our 
results with the parametrization of the data at a low normalization point performed 
recently by Gliick, Reya et al. [jl5|, 



few preliminary numerical results have been announced in ref. IT 
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2 From QCD to the effective chiral theory 



2.1 Light-cone representation for distribution functions 

The unpolarized quark distribution function of flavour / inside a nucleon with 4- 
momentum P, averaged over its spin, is given by the following QCD equation (see 

e.g. my- 



oo 



•(P|V;/(0)7+Pexp|-i5 dz'^A^{z')j^f{z)\P) (2.1) 
where '0/ are quark fields and is the gluon field. The antiquark distribution is 



oo 

1 



qj{x,^s) = -— J dz-e-'^P-^ 

— oo 

■{P\^f{Oh+Pexp{-ig rdz'"A^{z')]i;f{z)\P) . (2.2) 

L JO J z+=0, z±=0, fj. 

Here we use the light-cone coordinates 

and the nucleon state is normalized by 

{P\P') = 2P^6^{P -P') (2.4) 

where P is the nucleon 4-momentum. Throughout the paper x is the Bjorken variable, 
X = —(f'/{2P ■ q), where q is the 4-momentum transfer to the nucleon. 

The physical meaning of these equations is clear: the parton model description of 
the nucleon is justified in the infinite momentum frame; by using a backward Lorentz 
transformation one can recover a nucleon at rest but the price is that one has to work 
with the quark correlation functions on the light cone. The transition from the infinite 
momentum frame to the light-cone correlation function over the nucleon at rest is 
helpful since in the large Nc limit the nucleon is a heavy particle, and it is convenient 
to work in the nucleon rest frame. 

The right hand sides of eqs. ( p.l| , |2.2[ ) depend on the QCD normalization scale 
parameter fi. Actually, the non-local light-cone products of fields in these equations 
should be considered as a formal compact notation for an expansion in a series of local 
operators, each of these operators being renormalized at the scale /i. In contrast to 
QCD, the effective chiral field theory is nonrenormalizable and contains an explicit 
ultraviolet cutoff A (not to be mixed up with the Aqcd parameter). In the instanton 
vacuum model this cutoff appears as the inverse average instanton size: A ~ ~ 
600 MeV Q. The results obtained below refer thus to the low QCD scale fi of the 
order of 600 MeV. 
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We stress that we are computing the leading-twist distributions at a low normaliza- 
tion point and not the structure functions, observable in principle at low q^. The former 
differ from the latter by higher- twist power corrections which are large at low q^. The 
distributions we are computing can be directly used as initial conditions to the stan- 
dard perturbative evolution to higher where the power corrections are suppressed 
so that the distributions become directly related to the observables. 

Eqs.( p.l|j2.2| ) allow one to introduce a single distribution function qf{x) defined in 
the interval [—1, 1] identifying it at negative x with the antiquark distribution: 

(2.5) 

With this definition of qf{x) eq. ( |2.l| ) is valid in the interval — 1 < x < 1. 

An alternative way to the structure functions is via their moments. For example, 
one defines the moments of the singlet structure function as 




± 

I (2.6) 

-1 / 



Introducing local quark and gluon twist-2 operators [19| 



1 - 

Q<h^ll■■■^l^ = - — - [(jj-f^"^ D^'^ . . . D^'^ip + permutations - traces ] , (2.7) 

jn-2 

QG,^^l...^,„ ^ rj,^ [G'Mi!^2?M2 _ _ _ £)Mn-i(^^Mn + pgrmutatious - traces ] , (2.8) 
ni 

one can express the moments of the singlet structure function as nucleon matrix ele- 
ments of local operators: 

M„ = '^M^'^{P\i;f^{v^D^y--'^f\P), (2.9) 
where is a light-like vector, 

Vf.v" = 0, Vf.P'' = Mn. (2.10) 



When one starts to work with the effective low-energy theory (1.1) all information 
about gluons is already lost. Therefore, in order to recover the gluon distribution one 
has to go one step back, before the derivation of eq. (|1.1| ) Q and to rewrite the gluon 



operators in terms of (possibly nonlocal) quark operators. Ref. [16| explains how it 
can be done, and certain gluon operators are expressed there through quark fields only. 
Applying the method of ref. |16] to the twist-2 gluon operators (2.8) one observes. 



however, that in the leading order in the instanton packing fraction, 

/ = -p'^ « 1, (2.11) 

these operators are zero. The reason is the 0(4) symmetry of the instanton field: after 
integration over instanton orientations one can build the /ii . . . tensor only out of 
the Kronecker symbols, but it is impossible to get it traceless. In order to obtain a 
non-zero result one has to go beyond the zero- mode approximation of ref. fl^ and/or 
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consider contributions of many instantons. Both ways would lead to extra powers 
of the packing fraction of instantons. Meanwhile, it is the smallness of this packing 
fraction which is used in the derivation of eq. (1.1). Eq. (1.1) contains actually two 



dimensional parameters: the constituent quark mass M and the ultraviolet cutoff A; 
algebraically 

^ ~ / « 1. (2.12) 



Hence, to be consistent with the effective chiral lagrangian ( |1 . 1[ ) , one has to neglect the 
gluon operators, and to replace the covariant derivatives in eqs. ( |2.7| , |2.9|) by ordinary 
ones. That is what we are going to do in the rest of the paper. It should be kept in 
mind, however, that one has to introduce a finite cutoff A in order to make the nucleon 
mass and some of the structure functions finite (fortunately, the potential divergences 
are but logarithmical) . In fact gluons are resident in the "formfactor" of constituent 
quarks whose size is ~ 1/A. At A — > oo the constituent quarks are point-like, and 
there are no gluons. For the actual finite value of A ~ 600 MeV gluons have to show 
up. Unfortunately, the precise form of the gluon distribution depends, in the language 
of the chiral model, on the details of the ultraviolet regularization, therefore we shall 
not attempt to determine the gluon distribution in this paper. 

In the "quarks-antiquarks only" approximation one is left with the following ex- 
pression for the quark (and antiquark) distribution function qf{x) at — 1 < x < 1: 

oo 

Ifi^) = ^ I dz-e^^^^^-{P\i^f{^)^^^f{z)\P) (2.13) 



;+=0, 2x=0 



According to (|2.3D one can write 



V;/(0)7+V/(^) = V'}(0)7°7+V/(^) = ^^}(0)(1 + lW)^f{z) . (2.14) 



1 

72^ 

It follows from the constraint z+ = in ( |2.13| ) that = —z^. Taking the nucleon at 
rest we have = AI^/V^, where Mjy is the nucleon mass. Therefore, the quark (and 
antiquark) distribution function ( ^.13| ) is 



9/(^) = ^ / ^^°e'*"(p|^;(o)(i + 7V)^/(^)l^')L3=_,o,,^=o- (^-^^^ 

— oo 

For the polarized distributions, Agj(x), one has to insert the 75 matrix in this 
equation (see, e.g. |2^). In this paper we limit ourselves to two fiavours, u and d, 
neglecting altogether the strange quarks, therefore for both polarized and unpolarized 
distributions two cases have to be distinguished: fiavour singlet and flavour nonsinglet 
or isovector. For isovector distributions one has to insert the matrix in eq. ( |2.15| ). 
All four distributions are considered in this paper. 

Before turning to the calculation of these functions in the chiral quark-soliton model 
we make a general remark about the quark distributions in the limit Nc ^ 00. 

In this limit the quark part of the nucleon momentum in the infinite momentum 
frame is distributed among 0{Nc) quarks and antiquarks so that each quark and anti- 
quark carries 0{1/Nc) fraction of the nucleon momentum. This means that the quark 
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distribution function is concentrated at x ~ l/Nc- Keeping in mind that the total 
number of quarks minus the number of antiquarks in the nucleon is we can write 
the baryon number sum rule: 



J2 / dxqf{x)=N,. (2.16) 

/ - 1 

Comparing this sum rule with the fact that quark distributions are concentrated at 
X ~ 1/iVc we conclude that the singlet quark distribution has the following form in the 
large Nc limit: 

J2qf{x)=N^piN,x) (2.17) 
/ 

where the function p{y) is stable in the limit Nc ^ oo. Similarly, the isovector polarized 
distribution, Au{x) — Ad{x), is normalized, via the Bjorken sum rule, to gA which, 
theoretically, is of the order of Nc- Therefore, this distribution has also the form of 
eq. ( 2.17] ). It should be contrasted to the isovector unpolarized distribution, u{x) — d{x), 



normalized to the isospin which is 0{Nc)- Therefore, it is Nc times smaller than in 
eq. ( 2.17| ), and so is the singlet polarized distribution. In all four cases the antiquarks 



follow the same pattern as they are given by the same function as quarks but at x < 0. 
Exactly this behaviour of the distribution functions is obtained below. 

2.2 Chiral quark-soliton model of the nucleon 

Integrating out quarks in (^]^) one finds the effective chiral action, 

SesKix)] = -NcSplogD{U) , D{U) = ido - H{U), (2.18) 
where H{U) is the one-particle Dirac hamiltonian, 

H{U) = -i-i'^-i^dk + M-f°U^^ , (2.19) 

and Sp . . . is a functional trace. 

For a given time-independent pion field U = exp(i7r'^(x)r'') one can find the spec- 
trum of the Dirac hamiltonian, 

H^n = En^n. (2.20) 

It contains the upper and lower Dirac continua (distorted by the presence of the external 
pion field), and, in principle, may contain discrete bound-state level(s), if the pion field 
is strong enough. If the pion field has the unity winding number, there is exactly one 
bound-state level which travels all the way from the upper to the lower Dirac continuum 
as one increases the spatial size of the pion field from zero to infinity |^]. One has to 
fill in this level to get a non-zero baryon number state. Since the pion field is colour 
blind, one can put Nc quarks on that level in the antisymmetric state in colour. We 
denote the energy of the discrete level as -Eievj —M < Ei^^ < M. 

The limit of large Nc allows to use the mean-field approximation to find the nucleon 
mass - similarly to the Thomas— Fermi model of large Z atoms. To get the nucleon mass 
one has to add up NcE\q^ and the energy of the pion field. Since the effective chiral 
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lagrangian is given by the determinant ( 2.18| ) the energy of the pion field coincides 



exactly with the aggregate energy of the lower Dirac continuum, the free continuum 
subtracted. The self-consistent pion field is thus found from the minimization of the 
functional [Hi 



Mn = min N, I E,,, [U] + ^ {E^ [U] - ) L (2.21) 
[ £;„<o J 

From symmetry considerations one looks for the minimum in a hedgehog ansatz: 

C/c(x) =exp[i7r"(x)r"] =exp[i(n-r)P(r)], r = |x|, n=^ (2.22) 

where P(r) is called the profile of the soliton. 

The minimum ( |2.21 ) is degenerate in respect to the translation of the soliton in 



space and to the rotation of the soliton field in ordinary and isospin spaces. For the 



hedgehog field ( 2.22 ) the two rotations are equivalent. Including slow rotations of 
the saddle-point pion field and quantizing it gives rise to the quantum numbers of 
the nucleon: its spin and isospin components 0, |6|. In order to take into account 
the translational and rotational zero modes one has to make a unitary rotation of the 
quark eigenfunctions and to shift their centre, 



$„(x) ^P$„(x-X), (2.23) 

and to make a projection to a concrete nucleon state under consideration. The projec- 
tion into a nucleon state with given momentum P is obtained by integrating over all 
shifts X of the soliton, 

(P'l . . . I P) = 1 d^x e*{P'-P)-x , _ _ (2.24) 

The projection to a nucleon with given spin (5*3) and isospin (Ts) components is ob- 
tained by integrating over all spin-isospin rotations R of the soliton. 



(5 = T, S3, T3I . . . |5 = T, 53, r3) = j dR <t^V,T^{R) ■ . . 47iiR), (2.25) 

where </'f^(-R) is the rotational wave function of the nucleon given by the Wigner 
finite-rotation matrix |p: 

4;^^{R) = V2STT{-lf-'^-^D^_j,ls,^{R). (2.26) 

The four nucleon states have S = T = 1/2, with 5*3, = ±1/2. [Taking the next 
rotational excitation with S = T = 3/2 one can as well compute the A-resonance 
structure functions] . It is implied that dR in eq. ( ^^251 ) is the Haar measure normalized 



to unity, / dR = J d{AR) = J d{RB) = 1. In what follows we shall omit the superscript 
S = T. 
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2.3 Light-cone quark correlation functions in the nucleon 



In this subsection we derive several equivalent ways of presenting the quark distribution 
function ( p.l5| ) in the chiral quark-soliton model. Depending on the circumstances one 
can choose a more convenient representation. 

Eq. ( |2.15| ) is a matrix element of a non-local quark bilinear operator over the nucleon 
state with definite 4-momentum P and spin and isospin components. According to the 
previous subsection one can write down a general equation for such matrix elements; 
the time dependence of the quark operators is accounted for by the energy exponents. 
We write explicitly all the flavour {f,g = 1,2) and the Dirac = 1, ...,4) indices for 
clarity: 

(P,Ss,Ts\ip},{x^,^H''{y^,y)\F,S3,n) = 2PoN,l d^X J dR 4.^^-^^) 

■ J2 exp[i£;„(xO-yO)]<^,,(x-X)4^'i?^,ci>^'^(y-X)<A53T3(i?). (2.27) 

n 

occup. 

The functions are eigenstates of energy of the Dirac hamiltonian ( p.l9|) in the 
external (self-consistent) pion fleld Uc- Summation over colour indices is implied in the 
quark bilinears, hence the factor Nc in the r.h.s. 

In eq. ( ^.27 ) the quark is first annihilated in the nucleon by the operator V(y) 



and only then created by the operator '(/'^(x). Therefore the sum goes over occupied 
states. For the opposite ordering of the quark field operators, the order of creation and 
annihilation is opposite and the sum runs over non-occupied states: 



non— occup. 



We note that the eigenfunctions of the Dirac hamiltonian form a complete set of 
functions only when both occupied and non-occupied states are taken into account: 

^ <^lf^{^M^{y) = 5pl5{^ - y). (2.29) 

n 
all 

Adding up eqs. ( 2.27 , 2.28| ) at = and using the completeness condition ( 2.29| ) 
we observe that these equations are compatible with the standard equal-time anticom- 
mutator, 

{V},(x),^5^(y)} = 5f5f5(x-y). (2.30) 

In a Lorentz-invariant field theory the fermion operators should anticommute for 
any space-like separation: 

{^{y),^\x)] =0, if(x-y)2<0. (2.31) 

Our starting point, eq. ( p. IS]) , is a nucleon matrix element of a product of quark 
operators with a light-like separation which should be understood as a limit of space-like 



10 



separations. Indeed, the hght-like separation is obtained from the infinite momentum 
frame. As long as the momentum is large but finite one has a space-like separation. 
Therefore, the ip, ip^ operators in eq. ( p. 151 ) anticommute, and two alternative though 
equivalent representations ( p. 27 , 2.2^ ) can be written for the distribution functions; 
one sums over the occupied states, the other sums over non-occupied states of the 
Dirac hamiltonian. 

Finally, there exists another representation for the structure functions, this time 
through the imaginary part of the Feynman Green function in the background pion 
field. Let G(a;i,pi; W2)P2) be the Fourier transform of the Feynman two-point Green 
function 



G(.«.x;!,»,y) = -fe«,y|-^-i^|.«,x) 

= y| (.?, + MC/-^-) 3.^^,._,^;(^t,-,,)_,„ l^°. X) . (2.32) 
the free Green function being 

G('Hpr,P2) = (2vr)V^)fa-p,) ^,^+^^^^ . (2.33) 

In principle, one can expand the Green function ( |2.32| ) in powers of the pion field, 
U — 1, and its derivatives, dU. For the time- independent pion field U one has 

G{p1,Pi; p0,P2) = 27r<5(p?-p0)5(p?,pi,p2). (2.34) 

Using this Green function we derive in Appendix A the following representation for 
the singlet distribution function: 



J2qfix) = - Im I 27r6{p'+p'-xMM) Tr[5(/, p, p)(7° +7')]- (2-35) 

Since the Green function 5(p°,p, p) can be expanded in the pion field and/or its 
derivatives one can use this representation to make a quick estimate of the structure 
function, see section 7. 
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3 Singlet unpolarized distribution 



3.1 Sum over quark levels 

We start with the flavour singlet distribution function, u{x) + d{x), and the antiquark 
one, u{x) + d{x). The singlet case is simpler than the isovector one since one can ignore 
the rotation of the soliton. Indeed in that case the indices of the orientation matrices 
in eq. ( ^.27 ) are contracted, R^R = 1, and the integral over orientations R in eq. ( 2.271 ) 



becomes trivial and equal to unity. Taking the nucleon at rest, P = 0, Pq = Mn, and 
using eq. ( p. 271 ) for a specific light-cone separation of quarks operator as suggested by 
eq. (|2.15| ) we get: 



oo 



occup. 

•</.(-X)(l + 7V);-^{^(z-X)[3_^ (3.1) 



We remind the reader that the same equation gives, at x < 0, the antiquark distribu- 
tion, u{x) + d{x) = — {u{—x) + d{—x)). Passing to the momentum representation, 

^n(x) = j ^c^„(p)e'(P-), (3.2) 
and integrating over the coordinates of the nucleon center of inertia, X, we obtain: 



0<J 

^x)+d{x) = ^^^ I dz^e'^^'^'° y e-^^"^°(n|(l+7V)exp(-izV)|n) 
27r J ^ 

^ occup. 

:iV,MAr ^ (n|(l + 7V)<5(^n+p'-xM^)|n), (3.3) 



n 

occup. 



where |n) = <^„(p); taking the matrix element implies integration over 

However, eq. ( ^.3| ) can be understood in any other basis provided "p^ = —id/dx^ is the 

momentum operator. 



An alternative representation for the distribution function arises from eq. ( 2.2^ 
where summation goes over non- occupied states: 



7x(x) + d(x) = -^^^ /dzV^*^^^" y e-^^"^°(n|(l + 7V)exp(-izV)|ri) 
2tt J V 

^ non— occup. 

= -N,Mn y (n|(l + 7V)5(^n+p'-2;Mjv)|n), xE[-1,1]. (3.4) 



non— occup. 



Both representations, over occupied and non-occupied states of the Dirac hamil- 
tonian, should be absolutely equivalent were the theory finite. However, the effective 



chiral action (2.18) implies an ultraviolet cutoff: the nucleon mass itself and its struc- 



ture functions diverge logarithmically. In order to support the equivalence of two 
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representations, eq. ( |3.3| ) and eq. (p.4p, the ultraviolet regularization should be intro- 
duced in such a way as to preserve the anticommutativity of the quark fields at a 
space-like separation. In particular, it means that the completeness of the fermion 
states, leading to the equal-time anticommutator ( p. 29 ) must not be violated. For ex- 
ample, a relativistically-invariant regularization by the Pauli-Villars method preserves 



the equivalence of eqs. (|3.3| , |3.4| ) and other necessary properties of the distribution 
functions. On the contrary, a naive energy cutoff in the Dirac continua violates the 
completeness, and is hence unacceptable. 

In all above formulae we implied the vacuum subtraction: whenever we write 
(A^l^liV) with a local operator A, we mean in fact {N\A - (0|^|0)|iV). Without this 
subtraction (A^|^|A^) contains a divergence proportional to 5^^^(0). The vacuum sub- 
traction is effectively equivalent to the subtraction in eqs. ( p.3| , 13.4] ) of the appropriate 
sums over eigenstates \n^^^) with eigenenergies En'^ of the free Dirac operator Hq. 

For the free Dirac hamiltonian one has obviously En^ = iVIPnP + M'^. Therefore 
the delta function 6{xM]\f — En^ — Pn) vanishes if sign (x) = — signE'^'^^ Hence in 
eq. (|3.3| ) the vacuum subtraction term vanishes at x > (that is for quarks) whereas 
in eq. ( |3.4[ ) the vacuum subtraction is not needed at x < 0, that is for antiquarks. 

Both quark and antiquark distribution functions must be positive. At x > the 
function u{x) + d{x) is given by eq. (^) without vacuum subtraction. We can rewrite 
this equation as 



u{x) + d{x) 



NcMn f 



2lT 



y: [i>up)(i+7V)'j>n(p) 



p'^=xMM-En 



X > 0. 



occup. 



(3.5) 

This expression is explicitly positive because (1 -|- 7*^7^)72 is an orthogonal projector. 

Similarly, for the antiquark distribution one can write an explicitly positive expres- 
sion: 



u{x) + d{x) = — [u(— x) -|- d(— x)] 
N,Mn f d^ps_ 



27r 



jji^ [^l(p)(l+7V)^n(p) 



-xMn-Eu 



X > 0(3.6) 



non— occup. 



Although these formulae are explicitly positive one should keep in mind that they 
are ultraviolet divergent. In principle, one can regularize them so that the positivity 
is preserved. However it is not obvious a priori that one can find a regularization that 



simultaneously preserves both positivity and the equivalence of eq. (3.3) with eq. ( |3.4| ). 
We have observed that the Pauli-Villars method is favoured in this respect, too. 



3.2 Trace representation for distribution functions 

All the above representations for quark distribution functions are written as sums of 
diagonal matrix elements of certain operators over either occupied or non-occupied 
eigenstates of the one-particle Dirac hamiltonian H. Therefore one can easily rewrite 
these sums as operator traces. To this end we apply to eq. ( p.3| ) the spectral decom- 
position 
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n 

occup. 



\n){n\e-" '^-^ = j du 5{lo - H)e-"^ (3.7) 

and obtain 



u{x) + d{x) = N^Mn J diuSp[d{uj-H)6{uj+p^ -xMN){l + -f°-f^)]-iH ^Ho), 

— oo 

(3.8) 

where Sp . . . is the functional trace. 

Similarly, starting from the quark distribution function in the form of the sum over 
non-occupied states, eq. (|3.4|), we arrive to the following representation: 



u{x)+d{x) 



+00 

-NcMn j dw Sp 5{u: 

£;iov+o 



H)5{uj+p'' -xMN){l + rr) -{H^H( 



(3.9) 

Both formulae are valid at — 1 < x < 1. We remind the reader that, according to the 
results of the previous subsection, in eq. ( p.§| ) one has to make a vacuum subtraction 
at X < whereas in eq. ( |3.9| ) the vacuum subtraction is needed at x > 0. 

Symbolical as they may seem, eqs. (^^, 3^) give a practical way of computing the 
structure functions. To saturate the functional trace one can use any complete set 



of functions. For example, the quark eigenfunctions of the Dirac hamiltonian ( 2.20| ) 



may be used. In this basis the hamiltonian H is diagonal but the momentum p'^ is 
not. Another way is to use the eigenfunctions of the free hamiltonian, the so-called 



Kahana-Ripka basis |21|. The trace representation is also helpful in deriving general 
relations in a laconic form, see below. 



3.3 Baryon number sum rule 

Let us show that the baryon number sum rule is automatically satisfied in the above 
equations. To get the baryon number sum rule one has to integrate ( |3.8D over x from 
— 1 to 1. In the r.h.s. of that equation x enters through the product xMjv where 
Mjv = 0{Nc), so that in the large Nc limit we can replace this integral by the integral 
over the whole real axis of x, which leads to the following result: 



dx[u{x) + d{x)] = / duj Sp 6{uj - H){1 + j'^Y) - {H ^ Ho) . (3.10) 



Owing to the rotational hedgehog symmetry of the soliton the term j^j^ gives no 
contribution, and we are left with 

1 

J dx [n(x) + dix)] = N, Sp [e{-H + ^lev + 0) - 9i-Ho)] , (3.11) 
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which counts the number of the filled levels of the Dirac hamiltonian, the number of 
the levels in the free lower Dirac continuum subtracted. According to ref. it is the 
baryon number of the state. We have thus proved the sum rule 

1 1 

J dx [u{x) + d{x)] = j dx[u{x) + d{x) - u{x) - = N^B (3.12) 

-1 

where B is the baryon number of the state; B = \ for the nucleon. 

3.4 Moments of distribution functions 

We define the moments of the singlet structure function as 

1 

Mn= j (3.13) 

1 / 

Let us multiply (|3.8| ) by and integrate over x. In the large Nc limit one can 
extend the integration region to — oo < x < oo. Integrating the quark distribution in 
the form of eq. ( |3.8| ) we obtain the following representation for the moments: 



j dLoSp[5{uj-H){iu+p^)''-\l+j^-f^)\-{H^Ho). (3.14) 

-oo 

Note that we need the vacuum subtraction here since this result is derived by integrat- 
ing (3^) over both positive and negative x's. 

Similarly, the representation ( |3.9D based on the summation over non-occupied states 
leads to the alternative expression for the moments: 



l-n 



dujSp (^(w - i?) (a; +/)"-!(! +/^a) - (^H ^ Ho). (3.15) 



0^,3 ^, 



Both representations in fact follow from a third representation where the integration 
over uj goes along the imaginary axis, to the right of -Eiev Putting = iuj' we can 
write 



du)' 
27 



Sp 



uj' + iH 



(icu'+p3)«-i(l+/^3) 



(H^Ho). (3.16) 



Indeed, closing the lo' integration contour to the upper cut (corresponding to the lower 
Dirac continuum plus the discrete level) or to the lower cut and pole (corresponding to 
the upper Dirac continuum) one immediately reproduces eqs. ( 3.14 , p. 151) . Note that 
we assume that the uj' integration contour in eq. ( p. 161 ) is chosen so that the discrete 
level belongs to the occupied states. 

The fact that different deformations of the uj' integration contour in eq. ( 3.16| ) 
lead to the two representations in terms of summation over occupied and non-occupied 
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states respectively, gives another proof of the equivalence of these equations. We remind 
the reader that earlier we have derived this equivalence from causality, see subsection 
2.3. It illustrates a general connection between causality and analyticity. 

However, an arbitrary ultraviolet regularization of the theory may easily destroy 
both causality and analyticity. One has to check that a particular regularization does 
not induce new singularities in the uj' plane, and that it does not prevent one from 
closing the contours to the upper or to the lower cuts. For example, the popular 
proper time regularization would violate the last requirement. On the contrary, the 
Pauli-Villars regularization does not spoil the analytic properties, as well as casuality, 
see subsection 3.1. 



3.5 Momentum sum rule 

Let us derive the momentum sum rule for the quark distribution functions. It holds 
true only if equations of motion are satisfied, i.e. the pion field U which binds quarks 
in the nucleon is the minimum of the functional ( |2.21| ) which can be symbolically 
written as 

Mn = N, Sp [0(Siev -H + 0)F] -{H^ Ho). (3.17) 
The saddle-point equation for the self-consistent pion field reads: 



Sp[e(^iev + - H)6uH] = . (3.18) 

where djjH is an arbitrary variation of the Dirac hamiltonian under a variation of 
the chiral field U. Let us consider a particular (dilatational) variation of the chiral 
field, U{x) — > U[{1 + ^)x] with an infinitesimal ^ such that SU = ^x^dkU . Then the 
correspondent variation of the Dirac Hamiltonian is 

5uH = M^^^x^dkU^^' = ([x''dk,M^^U^'] = C ([x^dk,H] - i^^j^dk 



Inserting this variation into the saddle-point equation (3.18) and taking into account 
that 

Sp(^(^lev + - H)[x''dk, H]) = Sp([F, 0(Eiev + - H)]x''dk) = 0, 
we get a useful identity: 

Sp(0(i?w + - H)j''jkdk) = . (3.19) 

Owing to the hedgehog symmetry of the self-consistent pion field one can write a 
general tensor 

Sp(^(Eiev + - F)7%fc5z) = 6kiA (3.20) 



with the zero constant A because of eq. ( 3.19 ). 

Having made the necessary preparations, we turn to the second moment of the 
structure function M2. Using the representation ( |3.14| ) for the moments we can write 
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Mo 



iVeM^^Sp [9{Ei,, -H + 0){H + p^){l + 7%3^ 
N.M^'Sp \9{Ei,, -H + 0){H + 



(if ^ i/o) 
-Ho). 



(3.21) 



At the last step we have omitted the terms which do not contribute to the trace owing 
to the hedgehog symmetry of the saddle point chiral field. The identity (3.19) allows 
us to ignore the p^j'^j^ piece. Comparing eq. ( 3.21| ) with the expression for the nucleon 
mass, eq. ( p. 17 ), we obtain finally: 



Mo 



dx x[u{x) + d{x) + u{x) + d{x)] 



(3.22) 



meaning that quarks and antiquarks do carry the total momentum of the nucleon. 
The Pauli-Villars regularization is again privileged in that it does not destroy the 
momentum sum rule. 
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4 Isovector unpolarized distribution 



4.1 Soliton rotation 

It is easy to see that the isovector quark distribution function, u{x) — d{x), vanishes in 
the leading order of the 1/Nc expansion. In order to obtain a nonvanishing result one 
has to consider rotational corrections to the classical soliton. In the leading order of 
the 1/Nc expansion the rotation of the soliton is taken into account by the rotational 
wave functions (/ys^T^iR) ~ see eq. ( p. 26 ). In higher orders in 1/Nc one has to take into 



account that the functional integral (1.1) goes over a time-dependent chiral field, 

f/(t,x) = R{t)Uc{:x.)RHt) (4.1) 
For the rotating ansatz we have 

idt - H{U) = R [idt - H{Uc) + iR^R] R^ ■ (4.2) 



This leads to the following modification of the leading-order e qs. (|2.27 , 2.34| ), which 



takes into account the time dependence of the Green function ( 2.3^ ) 



(P = 0, 5 = T, S3, n\^l^\x)^iJiy)\P = 0,S = T, ^3, Ta) 
= 2MNi j d^X j dR(pl^j,^{R) 

•Tr {i?trfl(yO,y - X\[idt - H{Uc) + iR^R\-'\x^,^-X.)}^s,TAR), (4.3) 

where T is an arbitrary matrix in spin and isospin. Tr . . . denotes an ordinary trace in 
isospin and Lorentz indices, as contrasted to Sp . . . standing for functional traces; 

The integrand should be now expanded in the angular velocity R'^R which should 
be replaced by the spin operator S according to the following quantization rule [§, ^: 

R^R ^ ^S^t", (4.4) 



where I = 0{Nc) is the moment of inertia of the soliton, see below, eq. ( 4.16| ). At large 



Nc the moment of inertia is large, and the rotation may be considered as slow. 

4.2 Expressing the isovector distribution through a dou- 
ble sum over levels 

In the leading order in l/N^. one can simply neglect R^R in eq. ([4.3|). In this approx- 
imation one reproduces eq. (^.l|). However, neglecting the angular velocity altogether 
makes the nonsinglet distribution vanish. The first nonvanishing contribution to the 
isovector unpolarized quark distribution comes from the term linear in R^R in ( |4.3| ). 
Inserting this term into the general formula for the quark distribution function (|2.15| ) 
we obtain 

00 

u{x) - d{x) = '-^^ j Jd'xJdR 43^3 (R) 

^3 —CO 
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■Tj:{r^t^R{1 + 7°73)(20, z - X| 

— r|0,-X)| 

idt-H{Uc) J z^^=-z'^,z^=o 



1 



idt - H{Uc) 



iR^R) 



(4.5) 



Applying the quantization rule (|4.4| ) and introducing the orientation matrix in the 
adjoint representation, 



D,,{R) = lTr(i?tr-i?r^), 



(4.6) 



we find 



u{x) — d{x) 



Svr/ 



•Tr|(r''(l + 7V)(^°,z-X|— ^ 



dR4^T^{R)D3b{R)S^cl)s3T3{R) I d'X 
1 



H{U,) idt - H{U,) 



0,-X 



2:" =— 2:", z 



4.7) 



Let us first compute the rotational matrix element. Strictly speaking, this matrix 
element contains noncommuting operators D^i)[R) and 5", and one may worry about 
their ordering. However, due to the summing over nucleon spin the result does not 
depend on the order: 



E / dR<p\^T^{R)DMS''(I^S3T-AR) 



S3 



I dR4^T3(R)S''Dsb{R)ct>S3T3{R) = -h''\2T^^ 

S3 



(4.8) 



Therefore we get: 



u{x) - d{x) = -{2Ts)^^^ f dz'e'-^'^^' / d^x V Tr|r'^(l + 7°^ 

24-kI J J _ 



•(z°,z-X| 



idt - H{Uc) idt - H{U, 



^3 = _^0 z±=0 



(4.9) 



Since H{Uc) is time-independent (it coincides now with the hamiltonian H of the 
previous sections) we rewrite the matrix element as 



(zO, z - X|— ^ r'^— ^|0, -X) 
lOt — H lOt — H 



dw 
— ( 

2-n 



(z-X| 



X) 



(4.10) 



Noticing that (z — X| = (— X| exp[i(p • z)] and integrating first over X and then over 
we get finally 
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n(x)-d(x) = -(2T3)^^ / da; 



a=l 



uj-H uj-H 



(4.11) 



The result is of the same form as for the singlet structure function p.q ). However, 
in contrast to the singlet distribution which is a single sum over occupied (or non- 
occupied) levels (see eqs. ( ^.3| , ^^)), the isovector distribution is a double sum over 
levels. To see that explicitly we saturate the functional trace in eq. ( [4.11 ) by a complete 
set of functions, say, by the eigenfunctions of the Dirac hamiltonian ^n{v) iii the 
momentum representation. Then the integration over a; is performed with the help of 
the 5-function, and we get: 



u{x) — d{x) 



247r/ 



a=l m,n 



^>t(p)r«(l + 7V)^m(p) 



(27r)3 (p3 _ xMn - E^)ip3 - xMn - K) 

(4.12) 



The denominators here should be in fact understood with an ie prescription following 



from the original eq. (4.3): 



E„ 



Efn + iO for occupied states, 

Em — iO for non-occupied states. 



(4.13) 



We note that the wave functions (p) generically have singularities in the complex 
plane, therefore one cannot, generally speaking, close the contour of integration in p"^ in 
eq. ( 4.12| ) and reduce the integral to a residue of one of the denominators in eq. ( 4.12 ) 0. 



4.3 Isospin sum rule 

Let us now check the isospin sum rule, 



dx{u{x) — d{x) — [u{x) — d{x)]} = / dx[u{x) — d{x)] = 2T3. 



(4.14) 



As usually for general relations, the derivation of this sum rule is more laconic when 
one uses the symbolic expressions with functional traces. Indeed, integrating eq. ( 4.11| ) 
over X and replacing in the large Nc limit the integration limits [—1,1] by the whole 
real x axis we obtain: 



/ - = -(271,)§ / ^ Sp(r"(l + ,»,3)_2_ ^ ,4.15) 

1 a=l 



^In that way one would reco ver eq.(25) from a recent paper We conclude that their eq.(25) does not 
agree with our result, eq. ( 4.12 ). 
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The hedgehog symmetry allows us to drop the term 7*^7^. Taking into account that 
the nucleon moment of inertia is @ 



12 i 27r ^ uj- 

a=l 



H 



(4.16) 



we immediately reproduce the sum rule (414). One should keep in mind that both the 
moment of inertia and the isovector quark distribution contain logarithmic divergences, 
so that their regularization should be consistent with one another, if one wants to 
preserve the isospin sum rule. 



4.4 Gottfried sum rule 

The Gottfried sum rule follows from the assumption that the antiquark distribution 
in the nucleon is isotopically invariant: 



dx[u{x) — d{x)\ 



-1 



dx[u{x) — d{x)\ = . 



(4.17) 



This relation does not follow from any fundamental principle of QCD. Nevertheless, 
certain models assume this symmetry. Experimentally, eq. ( [4.17 ) is violated rather 
strongly |24]. Let us see what does the chiral theory predict for the r.h.s. of eq. ( [4.17| ). 
Integrating eq. (4.11) over x from —1 to we find: 



dx[u{x) - d{x)] = {2T3) 



247r/ 



^Sp[r'^(l + 7V)e(-u;-p')^^ 
27r L UJ — H 



OJ 



- Hi 
(4.18) 

It can be checked that this quantity vanishes in the leading order of the gradient 
expansion, therefore the Gottfried sum rule is satisfied only in the limit of very large 
solitons. For real solitons the expression ( 4.1^ ) is generally non-zero. 

The Gottfried sum rule has been analyzed in the context of the chiral quark-soliton 
model by Wakamatsu |25|. Several suggestions for the r.h.s. of the sum rule have 
been considered in that paper, however none of them coincides literally with the exact 
result ( p^ . 
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5 Isovector polarized distribution 

The polarized quark distribution function (see e.g. IQ) is given by 



oo 
— oo 

(5.1) 

where 5*3 is the spin of the nucleon in its rest frame. At negative x the function Aqf{x) 
has the meaning of the polarized antiquark distribution, 

Aqfix) = Aqji-x) (5.2) 

(note the opposite sign as compared to the relation ( |2.5D for the unpolarized distribu- 
tions!) 

One can easily check that in the leading order of the 1/Nc expansion only the 
isovector polarized distribution survives. We remind the reader that in the case of un- 
polarized distributions, on the contrary, the sm^/et distribution is large in Nc- Similarly 
to the derivation of eq. ( |3.3D from eq. ( |2.15| ) we can write 



oo 

Au{x) - Ad{x) = (2^3)^^^ / (iz°e"*^-^" J dR 4,s,iR) 

— oo 

• e-'''-''{n\R^T'R{l+j'j')j,eM-i^'p')\n)^ns,{R). (5.3) 

n 

occup. 

Computing the rotational matrix element with the rotational wave functions (piR) 
(see eq. ([4.8D ) we get: 



Atx(x)-Ad(x) = -(2r3)^^^ 

DTT 



• f dz^'"'^^'" e-^■^""°(n|T3(l+7V)75exp(-^zV)l^^)• (5-4) 



n 

occup. 



One has to take here = 1/2 for the distribution functions in a proton and = —1/2 
for those in a neutron. This sum over occupied states can be rewritten in terms of the 
functional trace: 



Au{x) - Ad{x) = --{2T3)NcMn 

dw Sp \6{uj - H)5{uj - xMn)t^{1 + 7°7^)75] - {H ^ Hq). (5.5) 
Integrating this equation over x we reproduce the Bjorken sum rule: 
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dx[Au{x) - Ad{x) + Au{x) - Ad{x)] = 2T3 qa, 



(5.6) 



where qa is the nucleon axial constant. In deriving this sum rule one has to keep in 
mind that the nucleon axial constant qa in the leading order in is given by the 
following functional trace p6|, p7[: 



9A 



^ / dcu Sp 5{uj - Hy-f'^-f^To -{H^ Ho) . 



(5.7) 



It is understood that one has to calculate this trace with finite pion mass and only 
then go to the chiral limit. 



6 Singlet polarized distribution 



For this structure function one should replace in eq. ( ^.3| ) by a unity flavour matrix, 
hence B)t^R is replaced by 1. For that reason the matrix element in eq. ( |5.3| ) is zero 
in the lowest order in the soliton rotation, and one has to expand the quark Green 
function to the first order in the angular velocity, R^R, as for the isovector unpolarized 
distribution. Combining thus eq. ( |5.3| ) and eq. (4^) we get 



Au{x) + Ad{x) = '^^^^' (2^3) j dz°e"*^'^^" y" d^x j dR^l^^^{R) 



•Tr((l + 7V)75(^°,z-X| 



1 



iR^R) 



dt-H 



|o,-x; 



23 = -20^ 2j^=0 



<PSsT-j.{R)- 



(6.1) 



Making again the quantization substitution R)R — > z5'"r°/(2/) and integrating 
over the soliton orientations R we get a representation similar to eq. (^ 



Au{x) + Ad{x) 



Sirl 



•IV((1 + 7V)75(^°,Z-X|— i 



idt -H idt-H 



|0,-X) 



;3 = _z0. 2 , =0 



(6.2) 



Performing the same steps as in subsection 4.2 we rewrite it through the functional 
trace: 



Au{x) +Ad{x) 



INcMn f dio 
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^Sp(^7V75<5(c^+/ 



xM. 



N) 



UJ 



H 



UJ 



(6.3) 



which is similar in spirit to eq. ( 4.11]) . Again, repeating the derivation of subsection 
4.2 we can write the singlet polarized distribution as a double sum over levels: 
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Au{x) + Ad{x) 



NcMni 
SttI 



E 



d^p H(p)(l+7V)75^m(p) 
(27r)3 (p3 _ ^Miv - - xMn - En) 



(6.4) 



understood with the same ie prescription as in eq. ( 4.13| ). 

Integrating eq. ( |6.3D over x we obtain the fraction of the nucleon spin carried by 
quarks: 



Au + Ad^l dx[Au{x) + Ad{x)] = ^ / ^Sp(7°7^75;^ r^-^) = , 

(6.5) 

which coincides, as it should, with the expression for the nucleon singlet axial constant 
g^^^ in the chiral quark-soliton model ||9|, |2^. The model calculation of this quantity 

gives g^^ ~ 0.36 the rest of the nucleon spin (at low g^!) is carried by the orbital 
moment of the constituents and of the distorted Dirac continuum Q . At higher values 
of an increasing portion of the nucleon spin is carried by gluons. 
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7 Expressing distributions directly through the 
sohton field 

All expressions for the quark distributions we have derived above are in fact certain 
functionals of the self-consistent pion field f7c(x) which binds the nucleon. It would be 
helpful to write down these functionals in a more explicit form. That can be done but 
at a price of additional approximations which, however, appear to be rather accurate 
in practice. 

A good place to start from is eq. ( |2.35 ) which relates structure functions to the 



imaginary part of the quark propagator in the background field of the soliton, C/c(x). 
Expanding the propagator (p. 321) in the derivatives of the pion field dUc up to the 



second order we obtain two contributions to the singlet distribution: 

, , N^MnM^ f £k f d^p 2tt5{p^ - xMn) 

l^Qfi^)- Im 



^ 2tt J (2^)37 (2^)4 (M2-(p-A;)2-iO)(M2-p2_ ^0)2 

•{(m2 -p2) Tr [i/-^5(-k)^c7~^^(k)(70 + 73)" 
+ TV [#i7-^5(_k)^[7-75(k)(/ + ^3)j|^ 

where we have introduced the Fourier transform of the pion field, 

U{k) = [ d\ e-^C''') [Uc{r) - 1]. (7.2) 



The integration over the quark loop momenta p can be easily performed: one first 
uses the ^- function to integrate over p^ , then one integrates over p^ taking the residues 
of one of the denominators. The condition that the poles in p^ lie on different sides of 
the integration axis so that one gets a non-zero imaginary part, is > xM]\j for x > 
and < xAIjy for x < 0. For x < we change the dumb variable k — k, so that 
the condition can be written in a common form as /c^ > |x|Mjv. We get 



E 



qf{x) « sign(x) ^^^^^' / (0 ^ {?7(k)[C7(k)]t} e{k^ - |x|M^) 
d2p^ M2 + , _ |x|M^(fc3 - |x|M^)k2 ^ 

(27r)2 (M2 + p^ + k2)2' " (A;3)2 > ^'■'^> 

The last integral is logarithmically divergent. We regularize it by the Pauli-Villars 
method: one has to compute the structure functions replacing the quark mass M 
Mpy where Mpy is the Pauli-Villars regulator mass, multiply the obtained result by 
M'^/MpY and subtract it from the original distribution. At the moment it is the only 
practical way of regularizing the structure functions we know of, which does not violate 
casuality and analyticity, see above. Note that in calculating static characteristics of 
the nucleon the requirements on the regularization method are not so restrictive. The 
value of the Pauli-Villars mass is fixed from the value of the Fj^ constant Q: 

F^ = 4N I iN / - In 



V (27r)4 (M2 + A:2)2 ^ Mly j {2ttY {M^y + k"^)-^ 47r2 M2 ' 

(7.4) 
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Integrating eq. ( |7.3| ) over p_L and performing the PauH-Villars regularization we 
obtain finally for the singlet structure function: 



sign(2;) 



In 



M2 + 



NcMnM^ 



PV 



M2) 



9{t' - \x\Mn) TV n7(k)[C/(k)] 



(M2^ + k2)(M2+k2) 



(7.5) 



We call it interpolation formula as eq. ( |7.5D becomes exact in three limiting cases: 
i) low momenta, \dU\ <C M, ii) high momenta, \dU\ ^ M, iii) any momenta but 
small pion fields, | log f7| ^ 1. Therefore, we expect that the interpolation formula has 
a good accuracy also in a general case. As compared to exact calculations involving 
summation over all levels the use of the interpolation formula gives an enormous gain in 
computing time: for a given profile of the pion field in the nucleon one has to compute 
numerically just three integrals. 

If the spatial size of the soliton is large, meaning that momenta k in eq. (|7.5| ) are 
small, one can neglect in the quark loop integral and get a simple formula: 



u{x) + d{x) ^ sign{x)F^MN J ^^^(^' " ^^^1) T> {u{k)[U(k)]^) (7.6) 

This is a remarkable equation. First, it shows that if the nucleon has the size ro 
the singlet structure function is concentrated around the values of x l/(rQMiy). 
Second, we notice that eq. ( |7.6| ) is non- analytic in the pion field momenta. It means 
that one cannot make a gradient expansion either for the structure functions or for 
their moments. This situation is new: all static characteristics of the nucleon admit 
expansion in the gradients of the pion field ||6| . By the way, it means that there is no 
reasonable way to extract structure functions from the Skyrme model (which contains 
just two and four derivatives of the pion field), even if one invents a way to identify 
leading-twist operators in that model. 

It is amusing that the second moment of the distribution ( |7.6D is 2/3 instead of 
unity. Indeed, integrating eq. ( [7. 6]) over x with x one gets for the second moment of 
the singlet structure function 

"'"1^ J (2^)3'^3Tr(c/(k)|t/(k)lt) (7.7) 

fc3>0 

where k'^ can be replaced by k2/3 for the spherically symmetric hedgehog field. At the 
same time the leading contribution to the large-size nucleon mass is given just by the 
kinetic-energy term of the chiral lagrangian. 



MN-^Jd^^Tr [5iC/(x)a,C/t(x)]. (7.^ 



Comparing eqs. ( |7.7| , 7.8 ) we see that in the large-size approximation the nucleon 
energy at rest computed from the total energy carried by its constituents is 2/3 that of 
its mass! Exactly this "2/3" paradox has been discovered some time ago in a similar 



situation by Lorentz |29[| who attempted to calculate the energy of the electron from a 
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charge distribution bound by some unknown forces of non-electromagnetic origin. He 
found an interesting relation: E = ^mc^. The paradox is resolved when one exploits 
the equation of motion for constituents. In our case eq. ( |7.8D has no minimum, however 
the full functional for the nucleon mass, including the discrete level, has. As shown 
in subsection 3.5 the use of the equation of motion for the pion field is essential in 
establishing the correct momentum sum rule. 

Eq. ( |7.6|) shows that in the leading order of the above expansion the singlet quark 
distribution coincides with the antiquark one. In order to obtain a nonvanishing result 
for the difference of singlet quark and antiquark distributions one has to go to the next 
(third) order in expanding the quark propagator. This time we perform the expansion 
in the coordinate space and get an elegant expression: 



- qf{x)] « ^ / d'y j^Ji e-^^'^-^J da J dp J d^ 
f °° 

•ei,-fcTV ([diU(y + <e3)][5,C/(y + PS^e;,)]'^ [duU {y + 7^63)] [C/(y)]^) . (7.9) 



It is remarkable that this result is consistent with the baryon number sum rule 
( |3.12| ). Indeed, for the solitons of a large size (for which the above expansion is justified) 
the baryon number coincides with the winding number of the pion field: 

B = ^J d^ye^JkTT [{U^diU){U^d,U){U^dkU)] . (7.10) 

Integrating eq. (^) over x one immediately obtains 

1 

/ dxJ2[qf{x)-qf{x)]=N,B (7.11) 
/ 



where B is given by the soliton winding number, eq. ( [/".lO ). 



Turning to the polarized distributions we remind the reader that in the leading order 
of the 1/Nc expansion only the isovector function Au{x) — Ad{x) survives. Similarly 
to eq. (|7.5|) we obtain: 



Au{x) - Ad{x) « 1 (2T3)^^£^^ / ^e{k' - |x|M^) Tr {mpm'r' 



m ■ 



M2 + k2 (M|,^ + k2)(M2 + k2) 

where is given in eq. (|7.3| ). For large-size solitons it can be simplified to 



(7.12) 



Au{x)-Ad{x)^\{2T^)F^MM J ^ ^(A;^ - M;v|x|)Tr (c/(k)[t/(k)] V) . (7.13) 



Note that eqs. ( [7.12| , 7.13| ) are even in x; the odd part of the Dirac continuum contri- 



bution to Au{x) — Ad{x) arises in the next order in the interpolation formula and is 
thus small. 
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Equations (7.5, |7.9|, 7.12) can be used to make a quick estimate of the structure 



functions, before one plugs into exact numerics which, in any case, is rather labori- 
ous. One should not forget, however, to add the contribution of the discrete level. 
To be consistent, a contribution of the discrete level with a Pauli-Villars mass to a 
distribution function f{x) should be subtracted: 

f^x) - /]r(x) - -^/]§;^(x). (7.14) 

8 Numerical results and discussion 

We have calculated numerically two structure functions surviving in the leading order 
in N(.: the singlet unpolarized distribution (section 3) and the isovector polarized one 
(section 5); in both cases we get quark and antiquark distributions separately. 

A variational estimate of the best profile of the pion- field soliton (see eq. ( p.22| )) 
has been performed in ref. P] yielding for M = 350 MeV 



r, 



P(r) = -2 arctg , tq ^ 1.0/M, Mn ^ 1170 MeV. (8.1) 

This profile function has a correct behaviour at small and large distances and is stable 
in respect to small perturbations. In our numerics we use the analytical profile ( |8.1| ): 
its difference with the exact solution is small. The value of the Pauli-Villars mass 
needed for the regularization is found from eq. ( |7.4| ) to be Mpy ~ 560 MeV. 

We estimate the Dirac continuum contribution using the interpolation formulae: 
eq. ( [7. 51 ) for the unpolarized singlet distribution and eq. ( 7.12| ) for the isovector po- 



larized distribution. In fact we have performed the exact computation of the Dirac 
continuum contribution by direct summation over all levels: as illustrated by Figs. 1,2, 
its deviation from the interpolation formula proves to be small. The description of 



exact computations will be a subject of a separate publication [30|. Since the reader 
may wish to repeat the calculations with his or her favourite set of parameters, we 
suggest the use of the interpolation formulae which give reliable approximations to the 
structure functions but can be computed in a few minutes on a PC. 

As stressed many times in the paper, the quark distribution q{x) is a single function 
defined for both positive and negative x's. At x < it gives actually the distribution of 
antiquarks, q{x) = —q{—x). The x-even combination, that is q{x) — q{x), x > 0, which 
is actually the baryon number density, is finite when one takes the ultraviolet cutoff to 
infinity. Moreover, it should not be regularized at all as it corresponds to the imaginary 
part of the effective chiral action (^]^). Probably, an ultraviolet regularization can 
be introduced which takes into account automatically the non-renormalization of the 
imaginary part of the chiral action. The original regularization by the momentum- 
dependent constituent quark mass M{p) 0] seems to satisfy this requirement |31]. In 



this paper, however, we mimic that regularization i) by the relativistic Pauli-Villars 
regularization for the logarithmically divergent x-odd part of q{x), that is for q{x) + 
q{x), ii) by making no regularization for the q{x) — q{x) distribution. In the case of 
isovector polarized distributions, the logarithmically divergent combination is [An(x) — 
Ad(x)] + [An(x) — A(i(x)], x > 0, corresponding to the g\ — gi structure function; the 
other combination originates from the imaginary part of the effective chiral action and 
should not be regularized. 
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Figures 1-5 present our results for the following distribution functions: 



Fig 


5-1 


x[u{x) + d{x) + u{x) + d{x)]/2] 


Fig 


5-2 


x[u{x) + — u{x) — d{x)]/2] 


Fig 


5-3 


x[u{x) + d{x)]/2; 


Fig 


;.4 


x[Au{x) - Ad{x) + Au{x) - Ad{x)]/2; 


Fig 


;.5 


x[Au{x) - Ad{x)]/2. 



In Figs. 1 and 4 we plot the contributions of the discrete level and that of the Dirac 
continuum (computed via the interpolation formula) as well as their sum separately. 
The exact calculation of the sum is also shown in Fig.l. It can be seen that the 
interpolation formula gives a good approximation to the exact calculations. The same 
is seen from Fig. 2. The discrete-level contributions are given in Appendix B. One can 
see that our quark and antiquark distribution functions are positive (Fig. 3), and that 
the contribution of the Dirac continuum is significant. 

It should be emphasized that the distribution of antiquarks arising from the discrete 
level (see eq. ( |B.5D ) is definitely negative (and sizeable!) ^ Positivity of the structure 
functions is restored only when one includes the contribution of the Dirac continuum. 
Indeed, the existence of a non-trivial discrete level is due to a strong mean pion field in 
the large-A^c nucleon. Taking into account only the discrete-level contribution to the 
structure functions means ignoring the degrees of freedom stored in the field creating 
that level, hence the "negative probability", q{x) < 0. Adding the Dirac continuum 
contribution makes the system complete, and all the probabilities become positive. We 
have demonstrated above that all the general sum rules are fulfilled only when one adds 
the discrete level and the Dirac continuum contributions together. 

There seems to be a lesson here for all variants of bag models. It is usually assumed 
that the three quarks in the bag give rise only to quark distributions, however they 
inevitably produce also negative antiquark distributions. The baryon number computed 
from the quark distributions only, turns out to be less than unity; to restore unity one 
has to subtract the negative distribution of the antiquarks. This circumstance is not 
often emphasized. To cure this decease one would need to add the structure function 
arising from the forces which keep the three quarks bound, in this case from the bag 
surface. 

To stress it once again: the discrete level produces not only valence quarks (in the 
sense used in the deep inelastic scattering phenomenology), whereas the Dirac contin- 
uum produces non-equal distributions of quarks and antiquarks (see Fig. 3), therefore 
its contribution should not be identified with the quark-antiquark sea of the DIS folk- 
lore. 

In the large Nc limit the nucleon is heavy, so we have neglected its recoil. For that 
reason the structure functions do not automatically go to zero at x = 1. However 
at X ^ 1/Nc the distributions behave as ~ exp(— const • Ncx); numerically, even for 
iVc = 3 all distributions computed are very small at x ~ 1. Nevertheless, we should 
caution those who may wish to reconstruct the structure functions from the moments: 
at n > A'^c the n-th moments become sensitive to the tail of the structure functions at 
X ~ 1 where 1/Nc corrections become 100% important, and the moments become thus 
absolutely unreliable. 

■^In the extreme case of a very strongly bound discrete level when it approaches the lower continuum, this 
level would not produce quarks at all - only antiquarks, but with a negative sign! 
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We remind the reader that our calculations refer to the leading-twist distribution 
functions at a normalization point of about 600 MeV. In order to make a meaningful 
comparison with the data one has to use our distributions as initial conditions for the 
standard perturbative evolution of the structure functions to higher values of where 
the actual data are available. This evolution takes into account the bremsstrahlung 
of gluons and also their conversion into quark-antiquark pairs. It is well known that 
the perturbative evolution makes the distributions more "soft" . Therefore, it can well 
wash out the quark excess at high x and cure the deficit at low x. This part of the 
investigation remains to be done. 

However, we can still make a comparison if not directly with the data, then with 
the parametrization of the data at a low normalization point ^ ~ 600 MeV performed 



recently by Gliick, Reya et al. [15, 17]. The standard perturbative evolution of their 



distributions describes well all the existing data at larger values of . It is known, 
though, that the evolution in the opposite direction is highly unstable. However, we 
believe that the low-point distributions suggested in refs. [15, ^ are reasonable, and 
we compare our curves to their in Figs. 1-4. 

Despite using the parametrization of the data at a low value of , the "experimen- 
tal" distributions appear to be more "soft" than the calculated ones. On the whole it 
looks as if we have determined the distributions at an even lower normalization point 
than that of refs. [^, 17 [. Since at these momenta are large it may require quite a 
short evolution range to move our distributions to those of [^, 17 [. 

The calculated isovector polarized distribution shown in Fig. 3 appears to be sys- 
tematically less than the parametrization of ref. [17]. We would like to make three 
comments on this discrepancy. First, we think that the parametrization of the po- 
larized distributions is less reliable than that of the unpolarized ones as it is based 
on less data with larger errors. In particular, the authors of [17| have assumed that 
Am(x) — A(i(x) = which is not confirmed in the model we are considering: this quan- 
tity appears to be of the same order as Ait(x) + A(i(x) of ref. [17[, see Fig.5. Second, by 
choosing the parameters of the soliton in eq. ( |8.lD we have, unfortunately, implanted 
a somewhat small value of the constant to which the distribution of Fig. 4 is nor- 
malized according to the Bjorken sum rule: it is qa ~ 0.96 instead of 1.25. Third, it is 
known [|9|, ^] that this channel is particularly sensitive to the 'i/Nc corrections which 
are altogether neglected in this paper. 

On the whole we get a reasonable description of several parton distributions without 
adjusting the parameters of the model to make a "best fit". 

Finally, we would like to comment on a recent work [22[ where quark distributions 
have been estimated in the Nambu-Jona-Lasinio model which, after certain simplifi- 
cations, is reduced to the chiral quark-soliton model of ref. [^ considered here. Only 
the contribution of the discrete level to the unpolarized structure functions has been 
studied in that work. As explained above, this approximation leads to a number of 
inconsistencies. For example, one obtains a wrong (negative) sign for the singlet an- 
tiquark distribution function, while the baryon number obtained from integrating the 
quark distribution only, is less than unity. As to the isovector distributions also con- 
sidered in [^, their basic eq. (25) does not seem to agree with our result given by 
eq. (FT^ ). 
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9 Conclusions 



At large number of colours the nucleon can be viewed as a heavy semiclassical body 
whose Nc "valence" quarks are bound by a self-consistent pion field. The energy of the 
pion field is given by the effective chiral lagrangian and coincides with the aggregate 
energy of the Dirac sea of quarks (the free continuum subtracted). Therefore, to 
compute the deep-inelastic structure functions in the large Nc limit, it is sufficient to 
calculate the quark and antiquark distributions arising from the discrete level occupied 
by quarks, and from the (distorted) negative-energy Dirac continuum. In contrast to 
all variants of the bag model, the completeness of the states involved guarantees the 
consistency of the calculations. Indeed, we have checked the validity of the baryon 
number, isospin, total momentum and of the Bjorken sum rules. We have also derived 
an expression for the r.h.s. of the Gottfried sum rule, which differs from the previously 
suggested ones. To our knowledge, it is the first time that the nucleon structure 
functions are theoretically calculated in a relativistic model which preserves all general 
properties of parton distributions. 

In the academic limit of a very weak mean pion field the Dirac continuum reduces to 
the free one (and should be subtracted to zero) while the discrete level joins the upper 
Dirac continuum. In that limit there are no antiquarks, while the distribution function 
of quarks becomes 5{x — 1/Nc). In reality there is a non-trivial mean pion field which i) 
creates a discrete level, ii) distorts the negative-energy Dirac continuum. As a result, 
the above 5-function is smeared significantly (but still in the range of the order of 
1/Nc), and non-zero antiquark distributions appear even at a low normalization point. 
It should be stressed that antiquarks come not only from the Dirac continuum but also 
from the discrete level, whereas the Dirac continuum produces nonequal distributions 
of quarks and antiquarks. 

As to the gluon distribution, it depends on the details of the ultraviolet regulariza- 
tion of the effective chiral theory (corresponding to the "formfactor" of the constituent 
quark), and for that reason we have not attempted to determine it here. Also, we have 
not tried to make a "best fit" to the parametrizations of the data |15, ^] by adjusting 
the parameters of the model, such as the constituent quark mass and the way one 
regularizes the effective chiral theory. In fact we have used practically the same set of 
parameters as deduced |Q, |6| from instantons. It is remarkable that we are able to re- 
produce the basic features of the distributions, despite neglecting completely the 1/A''c 
corrections. It would certainly be preferable to use our leading- twist distributions as 
initial conditions to the standard perturbative evolution to higher values of where a 
direct comparison to the data is possible, instead of comparing to the parametrizations 
of the data at a low normalization point, which to some extent is model-dependent. 

We have shown that, from the point of view of large all distributions can be 
divided into "large" and "small" ones. The large ones are the singlet unpolarized quark 
and antiquark distributions and the isovector polarized distributions. The isovector 
unpolarized and the singlet polarized are, parametrically, Nc times smaller, which 
seems to be confirmed experimentally despite that in reality Nc is only three. 

We have found that the structure functions are non-analytic in the pion field mo- 
menta, though one can still write exact expressions for the structure functions in the 
limit of large-size nucleons, see section 7. For arbitrary sizes we have derived interpo- 
lation formulae for the structure functions which allow one to compute the structure 
functions in a few minutes on a PC and reproduce the results of exact calculations to 
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good accuracy. 

We have restricted ourselves to the case of u and d quarks only, though the gen- 
eralization to three flavours is quite simple in the collective-quantization technique - 
see, e.g., ref. [^]. 

The methods developed in this paper can be easily generalized to higher-twist 
observables, like the (72, /it--- structure functions. They can be also used to estimate 
power corrections to the numerous structure functions at relatively low q'^. 
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A Expressing distributions through the Feyn- 
man propagator 

The Feynman Green function in a stationary background field is defined as 
Gi,-(x,y) =i(0|T{Vi(y)V5i(x)}|0) 

= ie(xO - y°) ^M-iEnix^ - y°)]$n,i(x)^nj(y) 

n 

non— occup. 
n 

occup. 

In the free case (no background field) it comes to 

The sign convention is such that the imaginary part of the Feynman propagator is 
positive. 

Let us rewrite the singlet structure function, say, in representation of eq. (0) 
through the Feynman propagator. To that end we divide the full range of integration 
in z'^ in eq. (|3.1| ) in two parts: from to +oo and from — oo to 0. Let us call these two 
contributions to the full structure function qi{x) and q2{x) respectively. In the q2{x) 
part let us change the integration variable — > —z^, so that the new variable also 
runs from to +oo. We have from eq. (^.1|): 

Ztt J Jq ^ 

occup. 

(A.3) 

= / d^X r e— V e*^"^''^„(-X)(70 + 7')$„(z0n3 - X) . 

ZTT J Jo n 

occup. 

(A.4) 

Comparing eq. ( [A.4| ) with the definition of the Feynman Green function ( [A.l| ) we 
see that it can be written as 



= J d'X e-«^^-^° Tr G{z'^^ - X, -z"; -X, 0)(7° + 7')- (A.5) 

The other part, eq. (A.3), does not directly fit into the definition of the Feynman 



Green function, however its complex conjugate does. Indeed, making the complex 
conjugation of eq. ( |A.3| ) and changing the integration variable X ^ X — z^rv^, we get 
exactly the r.h.s. of eq. (|A.4D , therefore q\{x) = q2{x). Using eq. ( |A.5|) we get for the 
full singlet structure function 
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Vg/(x)=-2 Im^^^ [d^X re-"*^^'^" Tr G(z°n3-X,-z0; -X, 0)(7° + 7')- 

^ ^TT J Jo 

(A.6) 

Let us now pass to the Fourier transform of the propagator. For the time-independ- 
ent background field one writes: 

= /^/^ '^'^^^^ -p^)5(rf,Pi,P2)e^(^^-)-^(^^-^). (A.7) 

Putting it in eq. ( [A.6D we integrate over X, which yields {2tt)^6{pi — P2). Finally, we 
integrate eq. ( |A.6| ) in z'^ and obtain 

V = 2 Im / A IV[S(/.P.p)(7°+7^)]^ 

f 27r i (27r)4i pO + ^3 _ j-jv/^ + ^ ^ 

We have made this derivation starting from eq. (|3.lD where the summation over 



occupied levels is used. We could as well start from the equivalent eq. ( p. 41 ) where 
summation goes over non-occupied levels. Repeating the same steps as above we 
arrive to eq. ( [A.8| ) but with the opposite overall sign and the opposite sign of iO in the 
denominator. Since the two formulas must be equivalent in a "good" renormalization 
scheme, it means that a non-zero imaginary part of the whole expression arises solely 
from the iiT6{p^ — xM^) pieces of the denominators in both cases. Thus, we obtain: 

Y^qjix) = - Im J 2^5(/+p3_xM^) TV p, p)(7°+7')], (A.9) 



/ 



which is eq. (|2!35| ) of the main text. 



B Bound-state level 

We present here the contributions of the discrete bound-state level to the singlet unpo- 
larized and to the isovector polarized structure functions. This is the cases where the 
discrete-level contribution is well-defined and in fact large. The other two structure 
functions considered in the paper are expressed through a double sum over levels, hence 
the contribution of the discrete level is not specific. 

The bound-state level occurs in the grand spin K = and parity 11 = -|- sector of 
the Dirac hamiltonian ( 2.19 ). In that sector the eigenvalue equation takes the form ^ 



/ McosP(r) - f - MsinP(r) 

1^ ^ - M sin P{r) -M cos P(r) 

We assume that the radial wave functions are normalized by the condition 




(B.l) 



^We change the sign of the 75 matrix and hence of the profile function P{r) as compared to ref. |^. The 
75 matrix is now that of Bjorken and Drell. The profile function is equal to — tt at the origin. 
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drr'[h'Q{r)+jt{r)] = l. 







We introduce the Fourier transforms of the radial wave functions, 



(B.2) 



oo oo 




where 



Rklir) 



-JiMk7 



-11 



2 r' fl d \ ' sin kr 



TT k^ \r dr 



(B.3) 



(B.4) 



The bound-state level contribution to the singlet unpolarized structure function 
can be simply obtained from the general eq. (p.3p. We get: 



[u{x) + d{x)Ui{x) =NcMn 



'^Ah\k)+3\k)-2 



xMn - Eie 

k 



-h{k)m 



|2:Mjv--Eicv| 



(B.5) 



Note that the r.h.s. is positive for all values of x, in particular at x < where eq. (|B.5| ) 
determines in fact the antiquark distribution. Since q{x) = —q{—x), it means that 



eq. ( B.5 ) gives a negative distribution of antiquarks at x > 0. At the same time it is 
easy to check by integrating eq. ( |B.5| ) that the baryon number sum rule is fully saturated 
by the discrete-level contribution only. It means that without the subtraction of the 
negative antiquark distribution the baryon number is not conserved. Simultaneously it 
means that the total baryon number of the Dirac continuum is zero, though locally in 
x the antiquark distribution from the Dirac continuum does not necessarily coincide 
with the quark one, see the dotted line in Fig. 2. 

The bound-state contribution to the polarized isovector distribution function in the 
proton is obtained from eq. 



lAu{x) - A<i(i-)1„, = jWoAfr 



\xMj^-Ei, 



k^ 



-2^^^^^i^^h{k)m 



(B.6) 
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Figure 1: The singlet unpolarized distribution, x[u{x) + d{x) + u{x) + d{x)]/2. Dashed 
hne: regularized contribution from the discrete level; dash-dotted line: contribution from 
the Dirac continuum according to the interpolation formula, eq. ( |7.5| ); solid line: the total 
distribution being the sum of the dashed and dash-dotted curves, dotted line: the exact total 
distribution; squares: the parametrization of ref. . 
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Figure 2: The baryon number distribution, x[u{x] 
bution from the unregularized discrete level, eq. 
contribution; squares: the parametrization of ref. [15]. 



+ d{x) — u{x) — d{x)]/2. Solid line: distri- 
).5|); dotted line: exact Dirac continuum 
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Figure 4: The isovector polarized distribution, x[Au{x) — Ad{x) + Au{x)—/S.d{x)\/2. Dashed 
hue: regularized contribution from the discrete level; solid line: the sum of the contribu- 
tions from the discrete level and from the continuum according to eq. ( |7.12| ); squares: the 



parametrization of ref. 
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Figure 5: The isovector polarized distribution of antiquarks, x[Au{x) — Ad{x)]/2. Ref. ||T7 
assumes this quantity to be zero. 
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